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Abstract 

Precise measurement of neutrino beam direction and intensity was achieved 
based on a new concept with modularized neutrino detectors. INGRID (In- 
teractive Neutrino GRID) is an on-axis near detector for the T2K long base- 
line neutrino oscillation experiment. INGRID consists of 16 identical mod- 
ules arranged in horizontal and vertical arrays around the beam center. The 
module has a sandwich structure of iron target plates and scintillator track- 
ers. INGRID directly monitors the muon neutrino beam profile center and 
intensity using the number of observed neutrino events in each module. The 
neutrino beam direction is measured with accuracy better than 0.4 mrad 
from the measured profile center. The normalized event rate is measured 
with 4% precision. 

Keywords: Neutrino oscillation, T2K, Neutrino beam, Neutrino detector, 
extruded scintillator, wave length shifting fiber 



1. Introduction 

To investigate the fiavor mixing of neutrinos [H [2] and the neutrino mass 
splittings, long baseline (over 100 km) neutrino oscillation experiments are 
currently running and also being prepared. In these experiments, it is ex- 
tremely important to measure the beam direction and intensity to ensure the 
stable neutrino production from a primary proton beam. 

T2K (Tokai-to-Kamioka) ^ is a long baseline neutrino oscillation exper- 
iment. An intense muon neutrino beam is produced by using the 30-GeV 
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proton synchrotron at J-PARC (Japan Proton Accelerator Research Com- 
plex) in Tokai. The proton beam impinges on a graphite target to produce 
charged pions, which are focused by three magnetic horns. The pions decay 
mainly into muon - muon-neutrino pairs during their passage through the 
96-meter decay volume. After traveling 295 km, the neutrinos are detected 
by the Super-Kamiokande (SK) detector [4J in the Kamioka Observatory 
(Fig. [T]). The goals of the T2K experiment are to measure oscillation param- 
eters with a precision of 5(sin^ 2^23) ^ 0.01 and S^Am^^) ^ 10~^ eV^ via z/^ 
disappearance and to investigate 9is through z/^ oscillation. 



Super-KAMIOKANDE 




Figure 1: Schematic view of the T2K experiment configuration. At the near detector 
location, there are two detectors: one is an off-axis near detector and the other is the 
INGRID on-axis near detector. 

T2K adopts an off- axis beam configuration [5] ; the beam center direction 
is 2.5 degrees away from the direction of SK so that the muon neutrino beam 
has a narrow energy peak at ^ 0.6 GeV, which maximizes the effect of the 
neutrino oscillation at SK and minimizes the background for the signal. The 
neutrino energy, however, varies as a function of the off-axis angle. Therefore, 
it is important to monitor and control the beam direction precisely; the 
beam direction is required to be controlled within ±1 mrad. In addition, 
monitoring of the beam intensity is important to ensure stable beam neutrino 
production. T2K has a mumon monitor downstream of the beam dump. 
It measures the beam direction and stability by detecting muons from pion 
decay for every bunch. Since the muon monitor detects only high energy 
muons which penetrate the beam dump, the phase space of parent pions 
covered by the muon monitor is much different from the one of pions which 
produce neutrinos to the near or far detectors. INGRID (Interactive Neutrino 
GRID) is an on-axis near detector and measures the neutrino beam direction 
and intensity by detecting neutrino interaction events. The covered phase 
space of the parent pions are much closer to the one for the off-axis neutrino 
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detectors than the muon monitor. Therefore, the measurement by INGRID 
is more directly connected with the T2K neutrino property. It was designed 
to provide daily measurements at the design beam intensity (750 kW primary 
proton beam power). 

This paper reports neutrino beam measurements using the INGRID de- 
tector based on the first two physics runs: Run 1 (Jan.-Jun. 2010) and Run 
2 (Nov. 2010 - Mar. 2011). During this time period, INGRID recorded more 
than 99.6% of delivered beam corresponding to 1.44 x 10^° protons on tar- 
get (POT). Using data from the same period, the T2K experiment observed 
indications of z/^ appearance [7]. 

In section 2 and 3, the design and basic performance of the INGRID 
detector are described, respectively. Details of Monte Carlo (MC) simula- 
tions and criteria for neutrino event selection are described in section 4 and 
5, respectively. In section 6, results of neutrino beam measurements with 
INGRID are summarized. 

2. Detector configuration 

INGRID is located 280 meters downstream of the meson production tar- 
get, where the spatial width (la) of the neutrino beam is about 5 meters. 
Therefore, it is designed to sample the beam in a transverse section of 10 m 
X 10 m with 14 identical modules arranged as two identical groups along 
the horizontal and vertical axis. Two separate modules are placed at off-axis 
positions ofl[ the main cross, as shown in Fig. [2] to monitor the asymmetry of 
the beam. 

Each of the modules consists of a sandwich structure of nine iron target 
plates and 11 tracking scintillator planes as shown in Fig. |3]left. They are 
surrounded by veto scintillator planes (Fig.[3]right) to reject charged particles 
coming from outside the modules. The dimensions of the iron target plates 
are 124 x 124 cm^ in the horizontal and vertical directions and 6.5 cm along 
the beam direction. The total iron mass serving as a neutrino interaction 
target is 7.1 tonnes per module. Neutrino interaction events are selected by 
reconstructing the track of charged particles generated by interactions in the 
iron target. The horizontal and vertical profiles are reconstructed from the 
number of observed events in each module. The beam center is measured 
as the center of profile. Each of the eleven tracking planes consists of 24 
scintillator bars in the horizontal direction glued to 24 perpendicular bars 
in the vertical direction with CEMEDINE® PM200. Each of veto planes 
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±5m 

Figure 2: INGRID on- axis near detector. The 16 identical modules sample the neutrino 
beam in a transverse section of 10m x 10m. The center of the cross, with two overlapping 
modules, corresponds to the designed neutrino beam center (6> = 0°). 




Figure 3: Structure of the module. It is a sandwich made of nine iron target plates and 
eleven scintillator trackers (left). The module is surrounded by the scintillator veto planes 
(right). 
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consists of 22 scintillator bars segmented along the beam direction. There 
are 9592 scintillator bars in total. No iron plate was placed between the 10th 
and 11th tracking planes. 

All the INGRID scintillator bars were produced at Fermilab [8j . The scin- 
tillator bars are made of polystyrene, infused with PPO (1%) and POPOP 
(0.03%), and are produced by extrusion in the shape of a rectangular cross 
section (1.0 cm x 5.0 cm). A white reflective coating composed of Ti02 
infused in polystyrene is co-extruded. One side of the rectangular face, far 
from the photo-detector, is painted with ELJEN® EJ-510. The length of the 
scintillator bars is 120.3 cm, 111.9 cm and 129.9 cm for tracking, top/bottom 
veto and right /left veto planes, respectively. A hole whose diameter is about 
3 mm at the center of the scintillator bar allows the insertion of a wavelength 
shifting (WLS) flber for light collection. 

The WLS flber, Y11(200)M by Kuraray \2\ is used for the light collec- 
tion. The diameter of the flbers is 1.0 mm and flts the active region of 
the photo-detector (1.3x1.3 mm^). The flbers are cut to the lengths of the 
scintillators and the cut surfaces are polished with diamond blades (Fiberfln 
Inc. FiberFin 4). One side of the flber is attached to a Multi Pixel Photon 
Counter (MPPC, Hamamatsu S10362-13-050C) \10\ [n\ with a specially de- 
veloped connector The other side is painted with EL JEN® EJ-510 to 
increase the light yield at the far side. 

The MPPC signal is transported to the front-end electronics by a co-axial 
cable (Hirose U.FL-2LP-068). In the electronics, the integrated charge and 
hit timing are digitized and recorded with a 2.5 photo-electron (PE) thresh- 
old. The bias voltage applied to the MPPCs is controlled by the front-end 
electronics with a precision of 0.02 V. A detailed description of the electronics 
can be found in fT3] . 

Figure [4] shows an example of neutrino event candidates in one of the 
modules. 

3. Basic performance of the detector 

INGRID identifies neutrino events by detecting tracks from muons. The 
hit efficiency for muon tracks is monitored with the beam induced muons 
incoming from outside of the detector. The mean light yield and timing 
resolution are monitored with cosmic-ray data. 
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Spill# 222B61, Time 12T12a931T 
Module B 



Side View 
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Top View 
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Figure 4: Typical neutrino interaction event candidate in one of the modules. A beam 
neutrino enters from the left. The size of the circles is proportional to the observed number 
of PE at scintillator bars, and black lines show the reconstructed tracks. 



3.1. Mean light yield 

The mean light yield per 1cm of a muon track is monitored with inter-spill 
cosmic-ray data for each scintillator bar. A typical light yield distribution 
of one channel is shown in Fig. [5| The distribution is consistent with the 
Landau distribution. The mean light yields of all the channels are shown in 
Fig. [gI The average mean light yield normalized to unit length is 24 PE/cm. 



3.2. Hit efficiency 

The hit efficiency is evaluated by muons with the following procedure; 
the muon track is reconstructed without using the hit information in the 
scintillator plane being evaluated, and then in that plane, channels expected 
to have hits from the track trajectory are checked whether they have a hit 
or not. 

Figure [7] shows the result of the hit efficiency measurement with the beam 
induced muons. Figure [8] shows the efficiency as a function of the track angle 
measured using cosmic-ray muons. The track angle is defined as the angle 
between the designed beam direction and the reconstructed track. The main 
reason for the inefficiency is because of the small gap between scintillator 
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Figure 5: Typical light yield distribution. Light yield is normalized by the path length. 



Mean light yield of each channel 
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Figure 6: Mean light yields. 
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bars, so the efficiency depends on track angle; a particle with small angle has 
more probability to go through the gap. As a result, the efficiency is smaller 
than that expected from PE statistics with measured light yield described in 
section 13. 1[ 




1.02 



Hit efficiency 

Figure 7: Hit efficiency for all channels. 



3.3. Hit timing resolution 

Hit timing resolution is estimated by measuring the time difference among 
hit channels for cosmic-ray tracks. Figure [9] shows the time differences of each 
hit channel from the average of all channels after the correction for differences 
in the readout cable length and the light propagation time through the fiber. 
The RMS is 0.9 nsec, which corresponds to the timing resolution if all the 
channels have the same resolution. The width of the primary proton beam 
bunch was about 30 nsec during RUN 1 and 2, so this resolution is sufficient 
for selecting the beam events. 

4. Monte Carlo simulations for the neutrino beam measurements. 

The Monte Carlo (MC) simulations for the neutrino beam measurements 
consist of three main parts. The first is a simulation of the neutrino beam 
production, which predicts the neutrino fiux and energy spectrum of each 
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Figure 8: Hit efficiency as a function of reconstructed track angle measured by cosmic-ray 
data. 




Figure 9: Time difference of hit channels from the average hit times for cosmic-ray tracks. 
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neutrino flavor. The second is a neutrino interaction simulation, which cal- 
culates the cross-section of the neutrino interaction and kinematics of flnal 
state particles taking into account the intranuclear interactions of hadrons. 
The third step is a detector response simulation to reproduce the scintillator 
light yield, and the response of WLS flbers and MPPC. 

^.7. Neutrino beam prediction 

To predict neutrino fluxes and energy spectra, a neutrino beam Monte 
Carlo simulation, called JNUBEAM [3j, was developed based on the GEANT3 
framework [M] . We compute the neutrino beam fluxes starting from models 
(FLUKA [13 US] and GCALOR [H]) and tuning them to experimental data 
(NA61/SHINE [ig and Eichten et al [19j). Energy spectra at the center 



and end of the horizontal modules are shown in Fig. [10} Because each module 
covers a different off- axis angle, the neutrino energy spectrum at each module 
location is different. The difference in the average neutrino energy between 
the center module and the end module is about 0.2 GeV. Energy spectra at 
10 m upstream from INGRID are predicted with the same procedure in order 
to simulate the background events from neutrino interactions in the wall of 
the experimental hall. 



xlO'' 



Neutino Flux 



center module 




1 



23456789 10 
Neutrino Energy [GeV] 



Figure 10: Neutrino energy spectrum predicted by JNUBEAM at the center and end 
modules. 
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4^2. Neutrino interaction simulation 

Neutrino interactions with iron in INGRID are simulated using the NEUT 
program hbraries [20j . The neutrino interactions in the scintiUator tracker are 
generated, taking into account the mass ratio. Combined with the neutrino 
flux predicted by JNUBEAM, we expect about 4 x 10'' neutrino interactions 
for 10^^ POT in the 14 modules consisting of horizontal and vertical rows. 
For the background event simulation, the neutrino interactions on an equal 
mix of carbon and hydrogen are generated uniformly in the upstream wall of 
the experimental hall. The wall mainly consists of sand, but for simplicity 
carbon and hydrogen are used in the simulation. The number of generated 
interactions is normalized by the beam induced muon events (described in 



Sec, 5. 1.5) in data. 



4.3. INGRID detector response simulation 

Detector response simulation was developed with the Geant4 framework [14] . 
The simulation includes a detailed geometry of the experimental hall. 

The energy deposit in each scintillator bar is simulated by Geant4 library 
and is converted to the number of PE at each MPPC. The conversion factor 
from the energy deposit to a number of PE is determined based on the 
measured light yields with cosmic-rays. The cross section of the scintillator 
bar is tuned to reproduce the hit inefliciency due to the dead region of the 
bar. The scintillator quenching effect is simulated using Birk's law with the 
value measured in [21j. Attenuation in the fiber is taken into account based 
on the measured attenuation length [22j. The response of MPPC, such as 
saturation due to the finite number of photo-diodes, is modeled based on test 
bench measurements [23j. 

The dimensions and mass of the iron target plates are implemented with 
the design value. For analysis, we make a correction for the measured mass 



difference for each iron plate, as described in Sec. 5.3 



5. Neutrino event selection 

The neutrino beam profile is reconstructed from the number of neutrino 
interaction events at each module. This section describes the selection pro- 
cedure for neutrino interactions and the systematic error. 
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5.1. Selection criteria 

A neutrino interaction event is identified by a long track from a charged 
particle generated by the neutrino interaction. First, pre-selections are ap- 
plied to reject accidental noise events. Then, tracks are reconstructed using 
hit information. After that, charged particles from outside of the module are 
rejected with the veto planes and the reconstructed event vertex is required 
to be inside the fiducial volume (FV). In these selections, each module is 
treated separately. The event selection criteria are described in the following 
sections. 

5.1.1. Event definition 

When there are four or more hits in a 100 nsec time window, all hits 
within ±50 nsec compose an event. 

5.1.2. Pre-selections 

A tracking plane with at least one hit in both x and y layers is defined 
as an "active" plane. We use a right hand coordinate in which z is along 
the beam direction and y is upward in the vertical direction. Events with 



three or more active planes are selected as shown in Fig. IT There is a 
discrepancy in the number of events with no active plane between data and 
MC. Figure [12] shows a MC event display of a typical event with no active 
plane. The first tracking plane is not counted as the active planes, but are 
used fro veto. Most of the events originate in low energy particles produced 
by the neutrino interaction in the wall of the experimental hall. We consider 
that the discrepancy is caused by uncertainties to simulate these low energy 
particles. 

After the selection with the number of active planes, light yields averaged 
over X or y layers of the active planes are required to be larger than 6.5 PE 
for both X layers and y layers as shown in FigfT3| inefficiency due to this 
selection is negligible for muon tracks. 

5.1.3. Tracking 

After the previous selections, tracks in x-z and y-z projection are re- 
constructed independently with a simple algorithm; first, hits in the most 
downstream active layer are adopted as the end point of the track. Then the 
track is extrapolated to upstream layers by checking the upstream hits. A hit 
is included in the track if the hit position is within two scintillator bars from 



the straight line extrapolated from the downstream hits. Figure 14 shows an 
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Figure 12: Examples of the MC event with no active plane. The red circle shows a hit by 
the particle and the blue circle shows the hit by MPPC noise. 
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Figure 13: Light yield averaged over active layers. Events with more than 6.5 PE are 
selected. 



example of a reconstructed track. Tracking efficiency is checked with cosmic- 
ray data and the efficiency is ^ 95% for cosmic-rays passing three scintillator 
planes. 
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Figure 14: Examples of the reconstructed tracks. The size of the circles shows the observed 
number of PE at scintillator bars, and black lines show the reconstructed tracks. 



The vertex position is reconstructed as the most upstream point of the 
track for each projection. The track angle is obtained by fitting the hits 



composing the track with a straight line. Figures [15] and [16] show differences 
between true and reconstructed x and z vertices, respectively, for MC events. 



The RMS for the x vertices is 2.7 cm. Figure [17] shows the distribution of 
3D angle between true and reconstructed muon tracks for MC events. The 
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RMS is 3.8 degrees. 

After the tracking, some badly fitted tracks are rejected by using the 
position difference of the vertex z between x-z and y-z projections. The 



difference is required to be within ±1 plane, as shown in Fig. 18 
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Figure 15: Differences between true and reconstructed vertex position in the x direction 
for MC events. 



5.1.4' Timing cut 

To reject off-timing events such as cosmic-ray events, only events of ±100 
nsec from the expected timing are selected (Fig, 19). The expected timing is 
evaluated with the primary proton beam timing [3j , the time of flight of the 
particles from the target to INGRID, and the delay of the electronics and 
cables. The event timing is deflned by the hit at the start point of the track. 



5.1.5. Veto and Fiducial Volume (FV) cuts 

Two selections are applied to reject incoming particles produced by neu- 
trino interactions in the upstream material, such as the wall of the experi- 
mental hall. First, events which have a hit in a veto plane or the flrst tracker 
plane at the upstream position extrapolated from the reconstructed track 
are rejected. Event displays of events rejected by the veto cut are shown in 



Fig. 20 After the veto cut, the flducial volume (FV) cut is applied. The 
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Figure 16: Differences between true and reconstructed vertex position in the z direction 
for MC events. 
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Figure 17: Angles between true and reconstructed tracks. 
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Figure 18: Difference of the z vertex position between x-z and y-z projections. 
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Figure 19: Time difference between measured event timing and expected neutrino event 
timing. Events within ±100 nsec are selected. 



21 



FV of each module is defined as a volume composed of the 3rd to 22nd of 
the 24 scintillator bars in the x and y directions, and from the second to the 
ninth tracker plane in the z direction. Events having a vertex inside the FV 



are selected as shown in Fig. [21] The events rejected by these selections are 
identified as 'beam induced muon' events. 



Figure [22] shows the vertex distributions in the x and y direction after all 
cuts. The fraction of the background in the selected events is 0.4%. 
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Figure 20: Event displays of rejected events by the veto cut. 
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Figure 21: Vertex x and y distributions. A volume composed of the 3rd to 22nd scintillator 
bars in the x and y directions is defined as Fiducial Volume (FV). 



5.1.6. Event selection summary 

The result of the event selection is summarized in Table [H The MC 
simulation includes neutrino interactions in the wall of the experimental hall. 
The number of neutrino interactions in the wall is normalized by the number 
of beam induced muons in the real data. The MC simulation reproduces the 
reduction in data well. 
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Figure 22: Vertex x and y distributions after all the event selections. 



Table 1: Event selection summary for 1.44 x 10^^ POT. MC sample is normalized by 
POT. The MC simulation includes neutrino interactions in the wall of the experimental 
hall. The number of neutrino interactions in the wall is normalized by the number of the 
bea m induced muons in the real data. 



selection 


Data 


(%) 


MC 


(%) 


1 


# of active planes > 2 


8.53 X 10^ 


100 


9.02 X 10^ 


100 


2 


PE / active layers > 6.5 


8.53 X lO*^ 


99.9 


9.02 X 10*5 


99.9 


3 


Tracking 


8.01 X 10*^ 


94 


8.40 X 10*5 


93 


4 


Track matching 


7.74 X 10^ 


91 


8.10 X 10*5 


90 


5 


Beam timing 


7.73 X 10^ 


91 


8.10 X 10*5 


90 


6 


veto cut 


3.30 X 10^ 


39 


3.30 X 10*5 


37 


7 


FY cut 


2.18 X 10*^ 


26 


2.17 X 10*5 


24 
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5.2. Selection efficiency 

The neutrino event selection efficiency as a function of true neutrino en- 



ergy is estimated by the MC simulation and is shown in Fig. 23, The ^20% 



inefficiency at high energy for the CC interaction is due to the events in 
which muons are produced with a rather large angle: for such events, the 
muon escapes from the module before it penetrates two iron plates. 
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Figure 23: Neutrino event selection efficiency as a function of true neutrino energy. 



5.3. Correction factors 

In order to compare the data with the MC expectation, some corrections 
are applied to the number of selected neutrino events. The correction factors 
are for the iron target mass and accidental MPPC noise. 

The neutrino interactions in the MC simulation are generated with the 
design mass of iron target plates, so a correction is applied to the neutrino 
interaction rate in each module individually to based on the difference of 
the design mass to the measured mass of the module: —1 ^ +1% for each 
module. 

Accidental MPPC noise hits sometimes results in mis-identification of 
the vertex. The dependence on the event selectionas a function of MPPC 
noise rate is estimated by a MC simulation where the PE distribution and 
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hit timing distribution are varied to reproduce data noise rates. According 
to this estimation and measured noise rate, the number of selected events is 
decreased by 3% with the existence of MPPC noise. The number of events 
in the MC simulation is corrected to account this effect. 

5.4' Systematic errors 

In this section the systematic error on the number of selected events is 
described. We describe both the errors related to the event selection criterion 
and the correction described in section 15.31 

The systematic error of the light yield selection is negligible because the 
threshold value for the selection (6.5 PE/active layer) is much smaller than 
the measured light yield from cosmic-ray. The systematic error of the track- 
ing efficiency is estimated by comparing the efficiency for several lengths 
of track between the data and the MC simulation. The difference is larger 
for the shorter track events and the maximum difference of 1.4% is taken 
as the systematic error. For the selection on the difference of vertex z be- 
tween x-z and y-z projections, the systematic error is estimated by looking 
at the change of the selection efficiency when varying the selection threshold. 
The difference (2.7%) between this and the nominal selection is taken as the 
systematic error. To estimate the systematic error associated with the FV 
selection, the number of neutrino events is re-counted as the fiducial volume 
is varied. In order to check the contamination of the beam induced muon 
from the edge, the fiducial volume is separated to center, middle and edge 
parts with the same volume. The difference of the event rate per unit volume 
is 1.1%, which is taken as the systematic error. The discrepancy between the 
hit efficiency measured with the beam induced muon event and that of the 
MC simulation is assigned as the uncertainty in the hit detection efficiency. 
The relation between the hit efficiency and the number of selected events is 
estimated by the MC simulation. A systematic error of 1.8% is assigned. 
The timing measurement is good enough that the systematic uncertainties 
resulting from the beam timing cut and out-of-beam events are negligible. 
There is a discrepancy between the beam induced muon event rate estimated 
by the MC simulation and that measured from the data. The change in the 
background contamination fraction from this discrepancy is estimated to be 
0.2%, which is taken as the systematic error for the beam-related background 
correction. 

The uncertainty of the iron mass measurement, 0.1%, is taken as the 
systematic error for the iron mass. The change in the selection efficiency due 
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to the time variation of the measured noise rate is estimated to be 0.7%, 
which is taken as the systematic error for the MPPC noise correction. 

Table [2] shows a summary of the systematic errors. The total systematic 
error is calculated as a quadratic sum of the errors and is 3.7%. In addition, 
there is 2% uncertainty for the normalization due to the systematic error on 
the POT measurement. 



Table 2: Detector systematic errors 



Item 


Error 


[%] 




Average light yield per active layer 


<0. 


1 


Tracking efficiency 


1.4 




Vertex z matching 


2.7 


FV cut 


1.1 




Hit efficiency 


1.8 




out-of-beam events 


<0. 


1 


Beam timing selection 


<0. 


1 


Beam-related background events 


0.2 




Iron mass 


0.1 




Accidental MPPC noise 


0.7 


Total 


3.7 



6. Results of the measurement 

6.1. Event rate and comparison with the MC expectation 



Figures [24] and [25] show the daily rates of the neutrino events and the 
beam induced muon events normalized by protons on target (POT). The 
beam induced muon events are defined as the events rejected by the veto cut 
or the FV cut in SecjS) They are measured typically with statistical errors 
of 1.7% and 1.1% each day, respectively. The average event rate is slightly 
decreased (-0.9%) on Feb. 4, 2011. This is considered to have been caused 
by a change in the primary beamline condition. The muon yield measured 
by the muon monitor showed a consistent tendency. The chi-squared values 
calculated from the average rates of the neutrino events (beam induced muon 
events) before and after Feb. 4, 2011 are 111.5 (122.1) for 125 data points 
and 26.4 (25.6) for 32 data points, respectively. The event rate in INGRID 
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before and after Feb. 4, 2011 remained stable within statistical error and the 
beam intensity is stable during each of these two periods. 

The data to MC ratio of the neutrino event rate is calculated to be 
1.06±0.001 (stat.)±0.04(detector syst.)±0.02(POT error). The uncertain- 
ties from the neutrino flux prediction and neutrino interaction model are not 
included in the systematic error. 
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Figure 24: Daily event rate of the neutrino events normalized by protons on target. 



Beam induced muon 
' Average 




X^/ndf 122.1/ 125 (-Feb. 4th 2011) 
X^/ndf 25.6/32 (Feb. ftth 2011~) 



May 



Jun. 



Nov.; Dec. 



2011 
Jan. 



Feb. May. 



Date 



Figure 25: Daily event rate of the beam induced muon events normalized by protons on 
target. 
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6.2. Beam direction 

The profile of the beam in the x and y directions is reconstructed with the 
number of neutrino events in seven horizontal and seven vertical modules, 
respectively. The observed profiles are fitted with a Gaussian function. The 
profile center is defined as the peak of the fit. Finally, the beam direction is 
reconstructed as the direction from the proton beam target position to the 
measured center at INGRID. 

In order to monitor the stability of the beam direction, the number of 
neutrino events is accumulated on a monthly basis. Figure 26 shows the 



observed profiles in April 2010. Black points in this figure show the number 
of neutrino events in each module and the red dashed lines show the fitted 
Gaussian. 



Black and blue points in Fig. [27|show the history of the beam centers in 
the X and y directions. All the points were stable well within 28 cm, which 
corresponds to the requirement of 1 mrad for the beam direction. Because 
beam direction was adjusted in November 2010, the beam centers in the y 
direction for later months are slightly shifted toward the center. 
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Figure 26: Neutrino beam profiles for x (left) and y (right) directions measured in April 
2010. 



The systematic error of the profile center measurement was estimated by 
a toy MC simulation. In the simulation, the number of events at each module 
is changed within the range of the total detector systematic error of 3.7%. 
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Figure 27: History of the neutrino beam centers. 



100,000 profiles are generated and RMSs of reconstructed center values are 
taken as the systematic errors; 9.2 cm and 10.4 cm for the x and y center, 
respectively. 

From the beam center measurement and the survey between the pro- 
ton target and the INGRID detectors, the average beam direction in x 
and y direction are measured as -0.014±0.025(stat.)±0.33(syst.) mrad and 
— 0.107±0.025(stat.)±0.37 (syst.) mrad, respectively. The beam direction is 
measured with a precision better than the requirement. 

7. Conclusion 

We have reported the muon neutrino beam measurement with the T2K 
on- axis near neutrino detector, INGRID, for the T2K Run 1 and Run 2 
data (1.44 x 10^° POT in total). INGRID consists of 14 identical modules 
arranged in a cross around the beam center. This configuration enables us 
to sample the beam in a sufficiently wide area to measure the beam center 
with a minimum of material. 

The neutrino event rate is measured on a daily basis and remains sta- 
ble within the statistical error, which is typically 1.7%. The data/MC 
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ratio of the event rate is measured to be 1.06±0.001(stat.)±0.04(detector 
syst.)±0.02(POT error). 

Beam centers in the x and y directions are measured on a monthly basis 
and are stable well within 28 cm, which corresponds to the required 1 mrad in 
angle. The angle between the expected and measured beam center direction 
in the x and y directions are -0.014± 0.025(stat.)±0.33(syst.) mrad and - 
0.107±0.025(stat.)±0.37(syst.) mrad, respectively. 

We conclude that the neutrino beam intensity and direction remain stable 
during the measurement based on the precise measurement with INGRID. 
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